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Abstract 



Graphene, the atomic-thin layer of carbon atoms, was first isolated on an insulating 
substrate in 2004 by two groups in Manchester University [U [5] and Columbia [3] • 
Those milestone experiments established the Dirac nature of the charge carriers 
in graphene. The same year, C.L. Kane and E.G. Mele predicted that intrinsic 
spin-orbit coupling in graphene, if strong enough, would lead to a novel state of 
electronic matter called the Quantum Spin Hall (QSH) state [HIS]- The QSH state 
is characterized by conducting gapless edge states circulating around an insulating 
bulk. Those edge states are protected from moderate disorder and interactions 
by a new topological invariant of the Z2 nature. While the strength of spin-orbit 
coupling is too weak in graphene, it was soon predicted and verified by transport 
experiments [3 [8] that the QSH state is realized in HgTe/CdTe quantum wells. 

In this manuscript, I will summarize some selected aspects of this huge field of 
research focused on Dirac matter including graphene and topological insulators. By 
Dirac matter, we have in mind various systems whose excitations obey a relativistic 
Dirac-like equation instead of the non relativistic Schrodinger equation. This report 
is mainly focused on the 2D topological insulators using graphene as a guideline. 

In chapter [T| the semimetallic character of graphene is derived and the symmetry 
protection of the Dirac points are discussed while chapters [2] and [3] are devoted to 
Chern insulators and QSH insulators respectively. 
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Chapter 1 

Dirac fermions in graphene 



Graphene, the atomic-thin layer of carbon atoms, was first isolated on an insulating 
substrate in 2004 by two groups in Manchester University [H [2] and Columbia 
[3]. Before those milestone experiments, it was already predicted that graphene 
should have a remarkable band structure [9] hosting Dirac fermions as low energy 
excitations \10\ lll j [T2]. Nevertheless it was commonly believed that a strictly (2D) 
two-dimensional carbon layer should also be unstable towards buckling or melting 
due to thermal fluctuations. Those experiments evidenced the existence of 2D Dirac 
fermions by the measurement of a very particular Quantum Hall effect [21 13] which 
is specific to relativistic carriers. It was also demonstrated that the density of such 
carriers can be tuned using a remote electrostatic gate, thereby realizing the first 
graphene-based field effect transistors [1]. Field effect transistors have a strong 
potential for applications in electronic devices, but they are also ideal systems to 
investigate the scattering properties of Dirac particles, including Klein tunneling 

[nmattSKis]. 

Here we would like to emphasize that the massless Dirac fermions in graphene 
are robust in many ways. Firstly they emerge at the level of non-interacting system, 
in contrast to other proposals of Dirac fermions at the nodes of d— wave supercon- 
ductors [17] or organic compounds [18] . Secondly, their zero mass character (gapless 
spectrum) is protected by the combination of time-reversal and inversion symme- 
tries. Thirdly, interactions are quite inefficient in opening a gap or disrupting those 
quasiparticles [HI [201 EI] • Finally spin-orbit coupling is also too weak to open a 
sizeable gap (see chapter [3| . 

In spite of this robustness, it is interesting, for both practical and fundamental 
purposes, to study how the Dirac points can be gapped out. In view of transistor 
applications, it is indeed mandatory to design narrow conducting channels by con- 
fining the carriers. Furthermore such conduction channels must be easily switchable 
between conducting (ON) and insulating (OFF) states. Hence it is crucial to un- 
derstand and control how to turn semimetallic graphene into an insulator. On the 
fundamental side, those issues are closely related to the problem of mass generation 
for relativistic Dirac quasiparticles. Because they carry an internal isospin degree of 
freedom, Dirac quasiparticles can gain a mass in different ways. The various masses 
are characterized by matrices in isospin space and can be classified according to the 
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symmetries they break or not. Increasing the number of internal degrees of freedom 
(that couple with the orbital motion) tends to create additional classes of insulators 
with very contrasted properties (22] . 

The physics of Dirac fermions in graphene has been extensively reviewed in far 
more details elsewhere |19| [20l I21j . The aim of this short chapter is simply to 
emphasize the concept of mass (or gap) generation in graphene [22] and to show 
that various insulating phases can be built (at least theoretically) from semimetallic 
graphene by adding proper perturbations. 

This chapter is organized as follows. We first show briefly how massless Dirac 
fermions appear as low-energy excitations of graphene. We then describe various 
microscopic models leading to finite gaps at the Dirac points and discuss the sym- 
metries of such external perturbations. We restrict this chapter to spinless Dirac 
fermions where only four types of masses are possible. This restriction is relevant 
either when the spin is frozen (by an in-plane magnetic field for instance) or when 
the spin is totally decoupled from the orbital motion (which is the case in real 
graphene owing to the weakness of the spin-orbit interaction [23^ I24j). More in- 
volved mass terms can appear when the spin-orbit coupling is included (chapter |3]) 
or when graphene is exposed to circularly polarized light. 



1.1 Massless Dirac fermions 

Here we briefly review the simplest tight-binding model for graphene (including 
only nearest-neighbor hopping amplitudes) , and derive the corresponding low-energy 
theory. We discuss the Dirac nature of the low-energy theory and the protection 
of the Dirac points by fundamental symmetries. At neutrality (when undoped), 
graphene is a 2D semimetal with two isolated Fermi points. Upon raising/lowering 
the Fermi level (by adding/removing electrons), a 2D metal with electron- like/hole- 
like carriers is generated, and a circular Fermi surface is formed. The Dirac nature 
of such carriers is revealed by their scattering properties at scalar potential steps, in 
particular at an ambipolar (pn) junction |25 tl26l[27] . This Klein tunneling has been 
reported experimentally [131 IHl 113 US] thereby providing a complementary proof 
of the Dirac nature of carrier beside the Quantum Hall measurements pi [3] . 



1.1.1 Tight-binding model of graphene 

Honeycomb lattice. Graphene consists of a honeycomb lattice of carbon atoms with 
two interpenetrating triangular sublattices, respectively denoted A and B (Fig. 



1.1.1). In this structure, each carbon atom has six electrons: two electrons fill- 
ing the inner shell Is, three electrons engaged in the 3 in-plane covalent bonds sp^, 
and a single electron occupying the pz orbital perpendicular to the plane. We are 
interested in the (2D )two-dimensional fluid formed by the pz electrons. Much of the 
physics of this 2D electronic system is described by the single orbital tight-binding 
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Figure 1.1: Tight-binding model for graphene. Red open (green filled) dots for A 
(B) sublattice. Orange thick arrows for Sa (a = 1,2,3). The distance between two 
sites is a = 0.142 nm and the surface of the unit cell is Aceii = 3-\/3a^/2. 



Hamiltonian: 

Ho = tY,c\ri)c{rj) = tY, Yl 4(r^ + (5a)c^(rA) + H.c, (1.1) 

a=l,2,3 

where the sum {i,j) runs over nearest-neighbor (NN) sites and r^, and t ~ —2.7 
eV is the hopping amplitude between the pz orbitals of two adjacent carbon atoms. 
The operator c(rj) destroys a fermion in the orbital pz at site r^, and is also denoted 
CA{ri) or CB(ri) depending whether this site belongs to A sublattice or B sublattice. 
The sum over runs over the A-sites which form a triangular Bravais lattice 
spanned by the basis vectors: 



ai = VSa Bx, a.2 = - i^VSer^ + 3ey j , 



(1.2) 



where a = 0.142 nm is the length of the carbon-carbon bond. The vectors da defined 
by 



^1, 



- (±V3ex + e 



y ; ' 



-a e 



:i.3) 



connect any A-site to its three B-type nearest neighbors (Fig. 1.1.1). The hopping 
matrix elements between next-nearest neighbors (and more distant atoms) are ne- 
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glected which is justified by the fact that those corrections are roughly ten times 
smaller than the main hopping t. 



Electronic band structure. Owing to translation invariance, the two-dimensional 
momentum k is a good quantum number. In order to diagonalize the Hamiltonian 



Eq. (1.1), we use the Fourier transformation: 



Cair^) = ^^e-^^-^^Caik), (1.4) 
where a = A, B is the sublattice index and is the total number of sites. After 



substitution of Eq.( 1.4), the Hamiltonian Eq. ( 1.1 ) becomes diagonal in momentum: 



^o = t J]7(k)4(k)cA(k) + H.c., (1.5) 



k 

where one has defined the quantity: 



^(k)= e^''-^'^ = 2e^^W2cos^^^ + e-^^-^ 



(1.6) 



a=l,2,3 

The electronic energy spectrum is given by 

^(k) = ±t|7(k)|, (1.7) 

which describes a valence band (minus sign) and a conduction band (plus sign) that 
are symmetric with respect to E = 0. Note that this zero of energy corresponds to 
the common energy of the atomic orbitals on sublattices A and B. The valence and 
conduction bands touch at isolated points of the first Brillouin zone (FBZ) obtained 
by solving the equation 7(k) = 0. Mathematically those touching points, also called 
Dirac points (for a reason explained in the next paragraph), span an infinite array 
of discrete locations of the reciprocal space. Of course there is a huge redundancy 
and all the points linked by a vector of the reciprocal lattice actually describe the 
same physical state. There are only two inequivalent Dirac points: 

k = ±K = ±^e,., (1.8) 

in the FBZ. Other solutions of the equation 7(k) = can be linked by a reciprocal 
lattice vector to one of those two solutions, and therefore describe the same physical 
state. 

Note that the existence of isolated solutions of 7(k) = 0, preventing the system to 
become gapped, is robust even if some crystal symmetries are lost and more hopping 
amplitudes are added. For instance additional second-neighbor hoppings will break 
the electron/hole symmetry discussed above, but will not affect the existence of Dirac 
points. Other perturbations, like a uniform anisotropic deformation on one type of 
bond, only shift the Dirac points and modify the conical dispersion around them 
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Figure 1.2: Electronic energy dispersion -E(k) = zbt|7(k)| of graphene for ky = 0. 



|18) . In fact the touching points are protected by more fundamental symmetries, 



namely inversion and time-reversal symmetries as we will show in detail later ( 1.1.3 ) 



Electronic filling. We now discuss how the electronic states are filled by electrons. 
Since the system has one electron per orbital (or site), the bands described above 
would be completely filled (both valence and conduction) if were ignoring spin. 
Due to the twofold spin degeneracy of each k-state, neutral graphene corresponds 
exactly to the situation of half- filling of the vr bands, and the Fermi level is at 
E = 0. Therefore the valence band is completely filled at zero temperature while 
the conduction band is completely empty, the Fermi level being reduced to the two 
contact points at k = itK. The total density of pz electrons, Uneutraiity = '^/^ceii — 
4.10^^ cm ~^ calculated as 2 electrons per unit cell (1 electron per carbon atom) is 
just the density necessary to fill completely the valence band. This density is not 
to be confused with the density of carriers n = utotai — nneutraiity which is zero at 
neutrality and can be positive or negative depending whether electrons are actually 
added or removed from the crystal. The carrier density is usually smaller than 
f^neutraiity, typically n ~ 10^^ — 10^^ cm~^ otherwise the Dirac Hamiltonian is no 
longer a good approximation for larger level of dopings (and for the corresponding 
high energies). 

1.1.2 Low energy theory near the Dirac points 

Dirac Hamiltonian. We consider now the low-energy theory for the single-particle 
states near the Dirac points. The momenta are written as k = itK + q close to 
the zero energy points (|q|a <^ 1), and the annihilation operators for those states 
are denoted ca±'k{^) = CA(iK + q) where q = qx^x + is a small momentum 
deviation from the Dirac points. From the first-order expansion of 7(k) around the 
Dirac points: 

7(±K + q) = -3o(±g^. + igy)/2, (1.9) 
5 



we find that the Hamiltonian describing the low energy excitations near k = 
= ±1) can be approximated as: 



where vp = —2>at/2 ~ 10^ m.s^^ ~ c/300 is the Fermi velocity. The Fermi velocity 
is basically the bandwidth t divided by the Brillouin zone (BZ) size 1/a. It is 
convenient to introduce a spinor representation: 

4(q) = (c!ik4kc!i-k4-k)' (1-11) 

embedding the two zero energy points = ±1. Then the single-electron Hamiltonian 
can be written in the compact form: 

4 

Ho = Y.Y. 4(q)(^o(q))af,C6(q), (1.12) 

q a,b=l 

with: 

UoicO = vriqxCrxTz + Qycry), (1-13) 

which has exactly the form of the Dirac Hamiltonian describing a spin one-half 
relativistic particles with zero mass. In particular the dispersion relation is simply: 

E{ci) = VFh\, (1.14) 

typical of a relativistic massless particle with velocity of light replaced hy vp. 

Nevertheless we would like to emphasize the differences between the Dirac equa- 
tion in the context of graphene system and in the high-energy framework respec- 
tively. In high-energy physics, the Dirac equation comes from Lorentz-invariance 
and very general considerations to associate special relativity and quantum me- 
chanics. Then the minimal objet to satisfy such an equation is a bispinor combining 
the spin and particle/hole symmetry relating positive and negative solutions of the 
Dirac equation. The coherent interpretation of the negative energy states led Dirac 
to the prediction of antiparticles which led to the discovery of the positron. 

In graphene, the origin of the Dirac physics is totally different. As we have seen, 
the spinors originate from a k.p expansion around special points of a particular band 
structure. Hence in graphene, there is no fundamental issue with the negative energy 
states that are just the valence band states (these states are in fact bounded from 
below by the bottom of the valence band) . Finally the emergent Lorentz invariance 



of Eq. (1.14) is only valid near the Dirac point, namely for wave vectors q located 
in a disk whose radius is far smaller than the inverse lattice spacing 1/a, whereas 
Lorentz invariance applies in the whole Minkowski space-time in particle physics. 

How many Dirac cones ? Fundamentally graphene is a two-band system because 
it has 2 orbitals per unit cell (one pz orbital per atom and 2 atoms in the unit cell), 
thereby having two states per momentum k in the first Brillouin zone (FBZ). Due 
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to fermion doubling, there are two species of massless Dirac fermions (one for each 
valley) carrying a sublattice isospin coupled to their momentum. Note that the 
presence of 4x4 matrices in Eq.(1.13) does not mean that graphene is a 4 band 
system in the same sense as the genuine 4 band insulators we shall study in chapter 
3} In deed for a given value of k in the FBZ graphene has only two states (Fig. 
1.1.1). Besides, when discussing transport or at least ballistic elastic scattering at 
the Fermi level, there are effectively 4 states sharing the same energy (Fig. 1.1.1). 
Then the physics depends on the ratio between intravalley and intervalley scattering 
rates. For instance, Klein tunneling and weak antilocalization are better observed 
if the intervalley coupling is much weaker than the intravalley coupling. 

Here we have not considered explicitly the real spin simply because it is not 
coupled to the momentum in the absence of spin-orbit. In fact at each valley, 
there are two completely degenerated and decoupled Dirac cones corresponding to 
each spin direction. Hence graphene has 4 Dirac cones with sublattice-momentum 
locking. This is at odds with surface state of 3D strong topological insulators which 
has a single Dirac cone with momentum coupled to the real spin. 



1.1.3 Microscopic symmetries 

In graphene, the fundamental symmetries are time-reversal symmetry T and inver- 
sion symmetry V. Again we restrict our discussion to spinless fermions. 

Time-reversal symmetry. In order to see how time-reversal symmetry acts in the 
effective low-energy theory, we first consider the wave function in the microscopic 
theory: 

*(r^) = ^^+(r)e^K-'--^ +(^^_(r)e-''K-'-^. (1.15) 

The time-reversed state is described by the complex conjugated wave function, 
namely: 

^*{va) = v\^{r)e-''^-'^ + ^\_{v)e''^-'\ (1.16) 

Now each Fourier mode UA±{<l)e^'^'^ of the envelope functions V3A±(r) is transformed 
into u^zp(q)*e~*^'''. Of course, one can proceed similarly for any site on the B sublat- 
tice. Hence in momentum representation, the effect of time-reversal transformation 
is to complex conjugate the amplitudes, change q into — q, and interchanges the 
valley labels while leaving unchanged the sublattice index. In the representation 
Ca(q) = (c^K'^BK'^yl-K'^B-K)' Operations described above can be summarized 
as: 

T = T,K„ (1.17) 

where Kc complex conjugate the amplitudes of the spinor and changes q into — q 
in momentum representation. It is very important to note that this time-reversal 
operation obeys = 1 because we deal with spinless fermions. When the spin is 
included (see chapter [s]) , the full time-reversal operation square to — 1 with Kramers 
degeneracy as a fundamental consequence. 
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Inversion symmetry. The inversion symmetry switches the sublattice A and B, 
and also changes the momentum k into — k. This means in particular that the 
valleys are switched and that the mode e***"" of the envelope function in valley ^ 
becomes the mode e~^'^^ in the opposite valley — ^. Hence the inversion operation 
V can be written as: 

V = a^T^, (1.18) 

and q — t- — q. 

Low energy theory Now we can check that the low-energy 4x4 Bloch Hamiltonian 
Eq.(1.13) is invariant both by the time-reversal operation 7~: 

TV.o{q)T'^ = vpTx [qxC^xTz + Qycrl) = ^o(-q), (1-19) 

and by the inversion operation: 

Vno{q)V'^ = vpCTxTx {qxCTxTz + QyCTy) a^T^ = V.o{-q), (1.20) 

where we have used the algebra of Pauli matrices. Using those symmetry operators, 
we can check the symmetry of any perturbation (or mass term) written in this low 
energy sector. 



1.2 Massive Dirac fermions 

We have seen that the combination of T and V protect the gapless nature of spin- 
less fermionic excitations in graphene. In 1984, Semenov discussed the effect of a 
sublattice staggered potential, which breaks the equivalence between A and B sites, 
thereby generating a gap at the Dirac points [11]. In a seminal paper [12j, Hal- 
dane introduced a special periodic magnetic field which also leads to massive Dirac 
fermions and further provides a lattice realization of the Quantum Hall Effect in the 
absence of a net magnetic flux through the unit cell. Finally it was recently empha- 
sized that a Kekule-type dimerization pattern of the carbon-carbon bonds can also 
open gaps at the Dirac points. In this report we call the corresponding insulating 
states the Semenov, the Haldane and the Kekule insulators respectively. In the low 
energy theory, a mass is a perturbation acting on the wave functions and able to 
open a gap. Such a perturbation is represented by a 4x4 matrix that anticommutes 
with the kinetic/ velocity part 'Ho{q) of the graphene Hamiltonian. 



1.2.1 Semenov insulator 

A local on-site staggered potential described by 

Hi = ^ eAC^Ai^A)cA{rA) + ^ £30^3(^3)03(^3), (1-21) 

spoils the equivalence between orbital energies ea and es on sites A and B re- 
spectively, and therefore also breaks the inversion symmetry V. This situation is 
naturally realized for a honeycomb structure where the A and B sites are actually 
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m_K = rriK 



I \ \ -k 

-K K 

tB{k) 




' \ \ -k 

-K K 

Figure 1.3: Low energy dispersion for the Semenov insulator and Haldane insulator. 



occupied by different atoms, like in BN crystals. Being local, this perturbation 
term is non dispersive (k-independent) and its Bloch Hamiltonian written is lattice 
isopin(8)valley space reads 

^i(k) =7^1 =eoI + Mia„ (1.22) 

which is typical of a massive relativistic particle. The mass is given by Mi = 
{ea — s_b)/2. The first term, proportional to identity I, only shifts the position of 
the Fermi level without changing the dispersion of the states. The o"^ = o"z (8) It 
term anticommutes with ?^o(q) and therefore can open a gap of size 2|Mi| at both 
Dirac points. The spectrum, obtained by squaring the Hamiltonian ?^o(q) + Afi^z 
and using {cr2,'Ho} = 0, is: 

E(q) = ±^vlq^ + Ml (1.23) 

There is a gap opening only if 7^ e_B, namely if A and B sites are non equivalents. 
We finally check explicitly that the perturbation 'Hi(q) is is odd under inversion 
'P = T^ax 

VHiV-^ = -Hi, (1.24) 

and even under time-reversal. 



1.2.2 Haldane insulator 

Haldane introduced a model of graphene under a modulated magnetic field that fully 
respects the spatial symmetries of the Bravais lattice but breaks T. This model was 
motivated by realizing the Quantum Hall effect without a global magnetic field by 
unit cell and without the Landau level structure. To realize that one should break 
time-reversal symmetry by inserting local fluxes which sum up to zero over each unit 
cell. This field preserves translational symmetry and the Bloch nature of electronic 
states. Those fluxes can be described by introducing unimodular phase factors in 
the second neighbor hopping amplitudes t2 — ?• t2e^^'^ ■, where the it sign corresponds 
to the different chiralities. 

The microscopic Haldane model is given by the Hamiltonian H = Hq + Hi + 



H2- The first perturbation Hi is the staggered on-site potential, Eq. (1.21), first 
considered by Semenov and discussed in the previous paragraph. Let us now focus 
here on the additional second-neighbor hopping term introduced by Haldane: 

H2 = t2Y,e"'^^'Pclcj, (1.25) 

where the sum runs over the next-nearest neighbor (NNN) sites {{i,j)). The chirality 
Uij = ±1 is defined as follows. Let us consider two NNN sites denoted i and j. We 

call d|j^ the unit vector from site i to the intermediate site (linking i and j), and 
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Figure 1.4: Tight-binding model for the Haldane model. Blue arrows (direction 
indicated by the oriented loops inside the hexagons) stands for t2e^'^. The vectors 
connecting NNN neighbor sites are defined as bi = ^2 — ^3, b2 = ^3 — <5i, and 
bs = 5i - S2. 



" (2) 

d - . the unit vector from this intermediate site to j. Then 



(1.26) 



To be more specific and to write explicitely the sublattice structure of those Haldane 
hopping terms, let us introduce the vectors bi = ^2 — ^3, b2 = ^3 — 5i, and bi = 
^2 — connecting next-nearest neighbor sites. Then the Hamiltonian H2 reads: 



H2=t2j2 ( I^c^(rA)cA(rA + b,)e''^ + ^4(rB)cB(rB + b,)e-"^ 1 + H.c, 

i V r.4 vb / 

(1.27) 

where i = 1,2,3. After Fourier transform, this Hamiltonian becomes 



n2(k) = 2t2 



cos{if) cos(k.bj)I -|- sin((/7) sin(k.bj)c72 

i=l,2,3 



i=l,2.3 



(1.28) 



which is valid for any k in the FEZ. 

The NNN perturbation is dispersive (k— dependent) because the it is nonlocal in 
real space. Near the Dirac points, one has simply to substitute k = K (or k = — K) 
as a zero order approximation. Here we focus on states in each valley and obtain: 



cos(K.bi) = sin(±K.bi) = 

1=1,2,3 j=l,2,3 



(1.29) 
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where we have used that K.bi = K.b2 = K.ba = 47r/3. The part of which 
is proportional to the identity just shifts the energies and spoils the electron-hole 
symmetry of the purely NN model. Nevertheless the system remains gapless (if 
ip = 0,tt) under introduction of a real NNN hoppings. This is because both V and 
7~ are preserved for real NNN hoppings. In contrast for complex hoppings, the term 
proportional to u^, denoted 'H2^{'k) hereafter, opens a gap at the Dirac points. In 
the low energy model, this Haldane mass term is jr2] : 

^^'^(±K) = -3\/3t2 sin(v?) a^r,, (1.30) 

which changes sign in different valleys. This is at odds with the Semenov insulator 
where both valleys are characterized by the same gap. The spectrum is again ob- 
tained by squaring the Bloch Hamiltonian ^^^{(l) and using the anti commutation 
property of the Pauli matrices: 



E = ±^Jvlp^ + 27t^sm'^ip. (1.31) 

This Bloch Hamiltonian breaks time-reversal symmetry because it is odd under 
the operator T: 

ml'^ici = o)T-^ = -nl'^{(i = 0), (1.32) 

but it is even under V = t^Ox- One can notice that this term cancels for = 
where the time-reversal is trivially restored. This is also the case at 93 = vr because 

„ — iTT _-| 



1.2.3 Kekule insulator 

The real- valued modulation of the nearest-neighbor hopping amplitude can open a 
gap when the wave vector of the modulation is connecting the two valleys. The 
Hamiltonian is + where the texture in the hopping amplitudes is described 
by [22]: 

iii = Y.Y. ^to.{rA)c\{vA^8o)cA{^A)^ii.c., (1.33) 

VA «=1,2,3 

where: 

<5t„(rA) = A(rA)e^K^"e^«-'-^ (1.34) 

The vector G = 2K connects the two Dirac points. The low-energy corresponds to 
the mass terms Tx^yCTz- 

The Semenov and the two Kekule masses belongs to the same triplet of the 
representation of Clifford algebra by 4x4 matrices while the Haldane mass is a 
singlet. This already signals that Haldane mass has a particular role that we shall 
explore in the next chapter on Chern insulators. 
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Conclusion and perspectives 



In this chapter, we have briefly reviewed that graphene has two Dirac points pro- 
tected by the simultaneous presence of time-reversal and inversion symmetries. 
Those Dirac points can be gapped out by perturbations that break one of those 
symmetries. We have discussed only spinless fermions in graphene but similar ques- 
tions can be raised for any semimetal where the carrier motion is coupled to one (or 
several) internal degrees of freedom. 

The number of distinct mass terms increases with the number of internal degrees 
of freedom that couple with the electronic motion. For the case of spinless fermions, 
there are four types of masses, namely the inversion breaking (Semenov) mass, the 
time-reversal breaking (Haldane) and two Kekule masses. Interestingly those insu- 
lating phases have very contrasted properties. The Semenov insulator is similar to 
an ordinary band insulator characterized by a purely local response to electric fields. 
In contrast the Haldane insulator has a chiral edge mode that conduct electricity 
non locally around an insulating bulk. The Kekule insulator can host fractional- 
ized excitations. Finally the inclusion of spin (and spin-orbit coupling) can lead to 
even more complex phases, including topological ones like the Quantum Spin Hall 
insulator (see chapter [3]). 

Finally, the general idea is that insulators are not all the same because they 
can differ by the topology of their Bloch wave functions even when their energy 
spectra are identical. This is at odd with usual non relativistic carriers in standard 
semiconducting systems where the gap is simple scalar and topologically protected 
surface states are absent. 
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Chapter 2 

Chern insulators 



Chern insulators (CIs), also called Quantum Anomalous Hall (QAH) phases, are 
band insulators that exhibit the Quantum Hall effect (QHE) in the absence of an 
overall external magnetic field. Historically the first CI model was devised by Hal- 
dane [12] who demonstrated that carefully chosen complex second neighbor hop- 
ping terms on the honeycomb lattice (chapter [T| generate a QHE carried by Bloch 
bands (instead of Landau levels). More recently the advent of topological insulators 
\28\ \29\ [30l l3T] has renewed the interest on the Haldane model and initiated the 
investigation of many models of CIs on various lattices [321 |33l [Ml EHl ES] . The 
Haldane insulator is also interesting in the doped regime where it becomes a metal 
with Berry phase effects and in particular where it exhibits the Anomalous Hall 
effect [37j. 

In order to generate a nonzero charge Hall conductance, it is necessary to break 
time- reversal (TR) symmetry. In the standard QHE, time-reversal symmetry is 
broken by a strong magnetic field which also creates the Landau level structure. 
In CI models the mandatory symmetry breaking is accomplished by introducing 
a periodic pattern of local fluxes that respects the full translational symmetry of 
the crystal. This makes CIs interesting in many respects. First, CIs differ from 
ordinary (atomic or covalent) band insulators by the presence of conducting edge 
channels. Second, CIs offer the opportunity to revisit the Quantum Hall physics 
in lattice systems which a priori differ from the usual continuous two-dimensional 
electron gas experimentally relevant for Si, GaAlAs and graphene-based quantum 
Hall devices. 

In practice generating Bloch bands with finite Chern number requires a very 
demanding band-structure engineering which is extremely difficult to achieve exper- 
imentally. In real materials, CI would require a nontrivial magnetic background, 
arising from a subtle interplay of spin-orbit (or some coupling between momentum 
and an internal degree of freedom) and exchange field effects (necessary to break 
TR symmetry). Another route could be to realize Chern (or QAH) phases under 
non equilibrium conditions [38l |39l SHI E] • 



This chapter is organized as follows. In Section 2.1 the general model of a two- 



band insulator with Berry phase effects in crystals is presented (42j . The section 



2.2 is devoted to the concept of (Chern) topological invariant and its relation to 
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the QHE in Chern insulators. Finally the mechanism for the formation of topo- 
logically protected edge states is explained by solving several simple models of an 
interface between two distinct gapped phases (section |2.3[ ). Conceptually CIs (or 
QAH phases) are the simplest topological phases, in terms of number of bands, that 
can be realized in a crystal: the simplest CIs are two-band insulators that provide 
building blocks for the understanding of more involved topological models like the 
QSH states (see chapter [s]). 



2.1 Two-band insulators 

Here we introduce the general framework to study Berry phase and topological 
effects in crystals using the simplest possible system: the two-band insulator (in 
mostly 2 dimensions). Two-band insulators play a fundamental role for the under- 
standing of the classification of electronic band structures, akin to the the two- level 
system in atomic physics. 



2.1.1 Bloch Hamiltonian 

Let us consider a d-dimensional lattice model for a two-band insulator. The tight- 
binding Hamiltonian in real space can be written as: 

H =Y1 4(ri)^ab(ri - rj)cb{rj), (2.1) 

where the sum over internal index a,b = 1,2 is implied. Ti'anslational invariance 
allows to use the d-dimensional momentum k as a good quantum number. Like 
for graphene (Eq 1.4), one introduces the Fourier transform of the second-quantized 
operators: 

ca{r^) = ^^e-^^■'^'ca{h), (2.2) 

where is the total number of lattice sites. This allows to write the Hamiltonian 
as: 

H = Y,ci{k)nab{k)cb{k), (2.3) 

k 

where ?^(k) is a 2 x 2 Hermitian matrix defined by: 

7^„b(k) = ^/i,6(r)e*^-^ (2.4) 

r 

This Bloch Hamiltonian is then a 2 x 2 Hermitian matrix, acting on the Bloch 
spinors, and that can be parametrized as 

?^(k) = eo(k)l2x2 + d(k).cr, (2.5) 

where cr = {cx, Cy, az) is the vector of standard Pauli matrices representing some 
internal isospin degree of freedom. This degree of freedom can be either a real spin, 
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the sublattice index (A and B sublattices of graphene), an orbital index (s and p 
orbitals defined on the same site). The details of the coupling is described by the 
vector d(k) = {dx(k.),dy {]<.), dz(k.)) of periodic functions of k. The structure of the 
Bloch Hamiltonian 'H(k) is constrained by the symmetries of the problem. 

All the information about the topology of wave functions is encoded in four real 
and periodic functions of the momentum, {eo{]s.),dx(k),dy{'k),dz{i^)), all defined 
on the whole Brillouin zone T^. The function eo(k) simply shifts the eigenvalues 
without affecting the eigenstates, and therefore it has no effect on the topological 
properties of the material. Nevertheless this function is very important eo(k) because 
it enters the spectrum dispersion and determination of the position of the Fermi 
level. 

Example: The Haldane model |12] is defined by the functions: 

3 

eo(k) = 2t2 cos((/?) cos(k.bi), (2.6) 

i=l 

3 

4(k) = me7(k) = t^cos(k.5„), (2.7) 

«=i 

3 

dy{k) = tlmj{k) = tY^sm(k.da), (2.8) 

a=l 

3 

4(k) = Ml + 2t2 sm{ip) ^ sin(k.bi), (2.9) 

1=1 

where the 2 dimensional momentum k lives in the first Brillouin zone T^. As 
seen in chapter [T} this model breaks simultaneously two symmetries, inversion and 
time-reversal. The terms eol2x2) CTxdxik), aydyCk) (which describes " unperturbated" 
graphene) are invariant under both T and V. In contrast the (iz(k)(j^ term is not 
invariant under T due to the flux term t2 siii{Lp) term [12], and not invariant under 
V due to the on-site staggered term Mi [TT] . 



2.1.2 Spectrum and wave functions 



The general Hamiltonian Eq.(2.5) can be easily diagonalized. The band structure 
consists of an upper (a = -|-) and a lower (a = — ) bands: 

i?„=±(k)=eo(k)±|d(k)|, (2.10) 

the corresponding wave functions being the spinors: 



a>+(k) 



u+i(k) 

U+2(k) 



cos f e-*^/2 
sin ^e^-^/^ 



(2.11) 
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in the upper band, and 



in the lower band. The k-dependent quantities ^ = and f = are the spherical 
coordinate angles of the unit vector: 



^cos (f sin 9^ 
nc^sin 
cos 6 



d(k)= I sinc^sin^ | , (2.13) 



which resides on the unit sphere 5^ while k spans the d-dimensional toroidal Bril- 
louin zone T*^. The mapping k — )• d(k) = d(k)/|d(k)| is essential and captures the 



topological properties of the Hamiltonian Eq.(2.5). We assume that the system is 
insulating so |d(k)| ^ everywhere and the mapping is well defined in the whole 
FBZ. 

2.1.3 Berry phase and anomalous velocity 

We can interpret k as a parameter that we can vary along a loop drawn in the 
FBZ and limit ourselves to d = 2. Along such a loop, the spinor will acquire a 
Berry phase which is the circulation of the Berry vector potential, also called Berry 
connection defined by 

2 

Aa(k) = i ^ Uaa'^kUaa, (2.14) 
a=l 

in each band a = ±1. The Berry curvature is the curl of the Berry connection: 

B,(k) = Vk A A,(k). (2.15) 
The fiux of Ba(k) through the whole FBZ (torus T^), 

Ca = ^ j dkBaik), (2.16) 

is called the Chern number of the band a. The Chern number is defined from the 
projectors on the lower (occupied) band. For the Haldane model, the valence and 
the conduction bands have a non zero Chern number (±1). 



Finally, from the Bloch Hamiltonian's expression ( |2.5[ ) one can define the com- 
ponents ji=x,y of the current operator: 
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which are 2x2 matrices for a two-band insulator. The second term, and in particular 
its non-diagonal elements, are responsible for Berry phase effects and eventually 
finite Hall response. 



2.2 Topological invariant 

The electromagnetic response of 2D two-band insulators can be computed using 



the Kubo formalism and the expressions of the anomalous currents Eq. (2.17). In 
d = 2 space dimensions, the Hall conductance is exactly given by the Chern number 
of the lower band. This is also the winding number of the mapping k — t- d(k) = 
d(k)/|d(k)|, which explains geometrically the quantization of the Hall conductance 
in such two-band model. This is reminiscent of the Thouless-Kohmoto-Nightingale- 
den Nijs (TKKN) invariant for quantum Hall systems |43j . 



2.2.1 Hall conductivity as a winding number 

The Hall conductivity can be calculated from Kubo formalism as [32]: 

where /±(k) are the occupation numbers of the conduction and valence bands. It 
is assumed that the Fermi level lies in the bulk gap. Hence at zero temperature, 
where /_ = 1 and /+ = 0, we have the relation: 

(^xy = -J^'^w^ (2-19) 

where is the winding number (or Pontryagin index) of the mapping k — )• d(k) = 
d(k)/|d(k)| between the Brillouin zone (torus T^) and the unit sphere (5*^): 

1 f .1, ( 5d(k) 5d(k)\ ^, 

= s / (i»t^ ^ j P-^"' 

In contrast to the Berry phase, this number is directly constructed from the pa- 



rameters d(k) of the Hamiltonian Eq. 2.5 (rather than from its eigenstates). This 



winding number is an integer that counts the number of times the unit vector d(k) 
wraps around the whole sphere S"^ while k is spanning the whole Brillouin zone 
In accordance with general classifications, there is a single number that character- 
izes the general structure of wave functions globally in k-space. This number is a 
relative integer, and it measures the charge Hall conductance in units of /h. To 
change it is necessary to change the parameter of the bulk Hamiltonian d(k) in 
such a way that the bulk gap closes. 

Interpretation as a magnetic textures We can see the field d(k) as some texture 
in momentum space. A terminology has been establish in real space to distinguish 
topological defects like skyrmions and merons. Those fields carry a topological 
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charge which is just the winding number defined above provided momentum k is 
replaced by position and also the FBZ is replaced by the manifold upon which the 
texture resides. 



2.2.2 Topological phase transition and band inversion 

Haldane model. Let us calculate this winding number for the simple model of massive 



Dirac fermions introduced previously. We use the parametrization Eq. (2.13) to 
rewrite the winding number as: 

1 r ine( de \ ^ 

^ At: J \dkx dky dk^ dky J 

= j t^^kV A (cos eVif) , (2.22) 

which in principle should be always zero because we integrate over a torus. Finite 
values can arise from singularities in Vy? that are always located at the poles 9 = 
and 6 = TT. Hence: 

n^ = -— / cosOV^p.dl, (2.23) 



Air „ 

where the integral is taken along loops encircling the poles. 

For graphene, the poles are reached when k is at the Dirac points. Then the 
sign of (iz(k = ±K) indicates whether the north or south pole has been reached. 
For the Semenoff mass, we have d^(k = ±K) = Mi) in both valleys. Then we have 
to notice that d = {^kx, ky, Mi) accumulates opposite phases while winding around 
the same pole (due to the presence of the valley index ^). For the Haldane mass, 
one has ^^(k = ±K) = MhC) which means that the accumulated phases (at south 
and north poles respectively) add up and finally: 

= (27r + 27r) sign{MH) = -sign{MH). (2.24) 
47r 

In fact for each valley, the winding of the angle ip is 2tt^ and cosO = Sign{M^) 
where is the mass in valley ^. Therefore one has 

nw = -^Yl '^<sign{Mi:) = ^ {sign{M.) - sign{M+)) , (2.25) 

which explains why the global winding number is zero when the masses are equal in 
both valleys, and why it is = ±1 in the Haldane phase characterized by a band 
inversion. This has been formulated in a more elegant and general way in Ref. 



Let us now consider the example of graphene in presence of some inversion break- 
ing and time-reversal breaking terms. So the mass matrix is (Mi — 3\/3t2 sin{ip)£^)az 
implying that the gap can close for Mi = ^3\/3t2 sinip in one valley (labelled by 
^ = ±1). This equality signals a one-electron topological quantum transition sepa- 
rating a QH insulator for a trivial atomic insulator. Finally we would like to make 
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Figure 2.1: Left panel: Interface between the Haldane insulator (x > 0) and the 
Semenov insulator (a; < 0). There is always a zero energy bound state located near 
X = (red dashed line). Right panel: dispersion E{ky) = sign{M H)fiVFky of the 
chiral edge mode. 



a comment on the terminology. This type of phase transition is purely a change 
between two one-electron Hamiltonian. It has in particular nothing to do with 
topological order defined by Wen. In particular the transition discussed here is 
not a transition between two topological orders. It is rather a transition between 
two band-insulators having distinct topological invariants (which characterize the 
winding of one-electron wave functions) . 

2.3 Chiral edge states 

Here we describe a rather general phenomenon occurring when an insulating 
(gapped) medium develops a region where the gap changes sign. Fermionic zero 
energy modes can show up in the gapless interfacial region, and get confined be- 
tween two fully gapped bulk regions. Two different situations are to be contrasted. 
First one might consider a mass kink without change of the topological invariant, 
for instance between two Semenov insulators, or two Haldane insulators (in the 
terminology introduced in chapter [2]) having opposite masses. Then the edge states 
exist but they are generally not protected. In contrast, the edge states residing at 
the interface between a Semenov phase and an Haldane phase (two topologically 
distinct insulators) is protected. We specialize our examples to 2D insulators so 
that the edge modes are running along ID interfaces. Nevertheless the idea is 
rather general and valid for D — 1 dimensional surface gapless modes emerging at 
interfaces between D dimensional gapped phases. This physics is reminiscent of the 
Jackiw-Rebbi model introduced in field theory [33] and of the physics of solitons in 
polyacetylene (D=l) [351136]. 
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2.3.1 Interface between topologically distinct insulators 

We assume that the half-plane 2; < is filled with a " Semenov" (inversion-breaking) 
insulator while an "Haldane" (time-reversal breaking) insulator occupies the half- 
plane X > 0. In principle one should define this heterojunction on the lattice by 
varying the parameters of the Haldane model (namely the on-site mass and the chiral 
phase (p) near the interface. Since we are mainly interested in eventual zero modes 
confined near the interface x = 0, we use the low energy effective model valid near 
the Dirac points (energies smaller than the bandwidth t). Using the translational 
invariance along the y— direction, the wave equation reads: 

{-ihvFCTx'rzdx + hvpkyay + M(x)) ^ = E^, (2.26) 

with M{x) = MsQ{—x)(Tz + MHQ{x)azTz- In fact we can even consider a more 
general shape by replacing the Heaviside functions Q{x) by smooth functions inter- 
polating between zero for negative arguments and unity for large positive arguments. 
Nevertheless the sharp step model is accurate provided the length scale for the vari- 
ation of the lattice parameters is smaller than the extension of the eventual edge 
state, namely hvp /max(Ms , Mh) ■ 

Let us first show there is always a solution at = and ky = hy solving the 
equation: 

hvpdx^ = -ia^TzMix)^, (2.27) 



obtained by multiplying both sides of Eq. (2.26) by iOxTz- For x > 0, we have: 



hvpdx^ = -GyMn^, (2.28) 

and the bounded solution (decaying at x — )• 00) is the eigenstate of ay with eigenvalue 
sign(MH) 

For X < 0, there is an additional valley matrix in the wave equation: 

hvpdx^ = -ayTzMs^, (2.29) 

and the corresponding bounded solution is the eigenstate of ay with eigenvalue - 
sign(^M5'). So the matching is possible, and there is a zero mode at the boundary, 
only if the two solutions above correspond to the same eigenvalue of ay, namely if 

sign{MH) = -sign{^Ms). (2.30) 

For any choice of the masses, this equality is always valid in one valley which is fixed 
by the relative signs of Mh and Ms- Therefore one obtains a zero mode which is 
polarized in the valley ^ = —sign{MsMj{). 

Now in order to obtain the wave function and dispersion E{ky) of this edge mode, 
let us restore finite energy E and parallel momentum ky in Eq.(2.26). Without 



further calculation one notice that the zero mode at fcy = is also eigenstate of 
hypkyay, and therefore its expression is still valid at finite energy and momentum 



21 




Figure 2.2: Left panel: Kink in the Haldane mass with the sign change occurring at 
X = 0. There is always twofold degenerated zero energy bound states located near 
X = (red dashed line). Right panel: dispersion E(ky) = sign{M u)hvFky of the two 
independent chiral edge modes (solid and dashed curves respectively). 



with the dispersion: 

E = —sign{S,Ms)hvFky = sign{MH)fiVFky. (2-31) 

The edge mode is chiral and shows up in the valley that is experiencing a mass 
inversion at the interface. In the limit of large Ms, the Semenov insulator can 
represent the vacuum. By reproducing this calculation for various orientation of 
the interface it is easy to demonstrate that the Haldane insulator is surrounded by 
a ID edge chiral edge mode that circulates clockwise if sign^MsMfj) is positive, 
and anti-clowise for negative MsMh- Note that if we assume that the vacuum 
is represented by a large positive Ms, then the sign of MsMh is simply the sign 
of Mff = — 3\/3t2 sin((^) which is set by the chirality of the flux pattern in the 
microscopic Haldane model (see chapter [2]). 

2.3.2 Kink in the Haldane mass 

As was the previous one, this section might be a bit academic since it requires to 
make a junction between two Haldane insulators with opposite chiralities while there 
is not yet any experimental evidence of a Haldane phase in graphene. Nevertheless 
we think one can learn a lot from those simple toy models. 
The full wave equation for the Haldane kink reads: 

{-ihvpfTxTzdx + hvpkyay + MH{x)azTz) ^ = E^!, (2.32) 

where Mh{x) is a real monotonic function describing a kink with Mh (oo) positive 
and Mh{—oo) negative hereafter (the opposite case can be treated similarly). We 
take the origin x = where M{x) has its zero. We expect that a bound state might 
show up near x = because the insulator becomes "locally" gapless there. 
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We first look for a E = solution at ky = by solving the equation: 

ihvpcT^T^dcc^ = MH{x)a^T^-^. (2.33) 
By multiplying each side by —iaxTz, it is obtained: 

hvpdx^ = -MHix)ay^, 



which has the solution: 



^'(x) = exp ^- j dx'MH{x) 



)/hVF I \(Ty = +1) 



exp ( — j dx Mh{x) /hvp 



i , 
1 



(2.34) 

(2.35) 
(2.36) 



Hence there is a twofold degenerate zero mode at ky = 0. 

Now we can restore a finite transverse momentum ky and observe that the above 
solution is an eigenmode of hvpCTy with energy E = hvpky. The two degenerate 



chiral zero modes Eq.(2.36) yield two degenerate chiral modes propagating in the 
same direction along y-axis. This is consistent with the fact that the Haldane 
model breaks time-reversal symmetry. We can understand the Haldane kink as two 
remote Haldane insulators (with opposite chiralities) that would have been brought 
in contact together adiabatically. After such a process one would have to modes 
running in the same direction along the interface considered. 

2.3.3 Kink in the Semenov mass 

One can easily reproduce the similar analysis for a kink of the Semenoff mass by 
solving the wave equation: 



{-ihvFCTxTzdx + hvpkyay + Ms{x)az) ^ = E'^, 



(2.37) 



where Ms{x) is a real function satisfying Ms(0) = 0, AIs{oo) > and Ms{—oo) < 0. 
The equation for the eventual zero energy mode at ky = is then: 



hvpdx^ = -Ms{x)ayTz'i!, 



whose solution reads: 



^'(a;) = exp ^- j dx Ms{x') /hvp ] WyT^ = +1) 
dx Ms{x') /hvp 



exp 



/ 


n\ 






\ 




i 









a 





+ b 


1 




\ 











(2.38) 

(2.39) 
(2.40) 
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Figure 2.3: Left panel: Kink in the Semenov mass M{x)az with the sign change 
occurring at x = 0. There is always twofold degenerated zero energy bound states 
located near a; = (red dashed line). Right panel: dispersion E{ky) = zizhvpky 
of the two counter-propagating edge modes (solid and dashed curves respectively). 
This counter propagation is the natural consequence of the time-reversal invariance 
of the system. 

As a major difference with the Haldane kink, the two parts of the wave function leads 
to opposite chiralities when a finite ky is restored. This is because they correspond 
to eigenmodes of dy with opposite eigenvalues ±1. This is consistent with the global 
time-reversal symmetry of the system. 

Conclusion 

We have taken advantage of the conceptual simplicity of CIs to introduce the general 
ideas topological invariants and topologically protected gapless edge modes using the 
Haldane model as a guiding example. At the interface between two insulators with 
opposite masses (same matrix but sign inversion through the interface), or between 
two topologically distinct insulators, zero modes generically show up and form an 
interfacial metal. When the two insulators have distinct topological invariants, the 
metallic (gapless) of this interfacial electronic system (holographic) is topologically 
protected . This is quite similar to the situation of the chiral edge states of the 
integer QHE protected by the finite TKKN invariant. The ideas discussed in this 
chapter will be used in the following to understand other topological states, like the 
Quantum Spin Hall (QSH) state (see Chapter [s]) and Floquet topological insulators 
(see chapter ??). 



24 



Chapter 3 

Topological insulators 



In the absence of spin-orbit coupling, the Dirac points of graphene are protected 
by the combination of two fundamental discrete symmetries: time-reversal and spa- 
cial inversion (chapter [T]). In 2005, C.L. Kane and E.G. Mele demonstrated that 
the situation is drastically changed when the spin is coupled to electronic motion: 
intrinsic spin-orbit coupling does open a gap at the Dirac points without breaking 
any of those fundamental symmetries [H U]. The resulting insulator, the so-called 
Quantum Spin Hall (QSH), is a novel state of electronic matter that cannot be adi- 
abatically connected to a trivial atomic insulator without closing (and re-opening) 
the bulk gap. The QSH state is distinguished from ordinary band insulators by the 
presence of a one-dimensional metal along its edge [H f5| which is topologically pro- 
tected by a Z2 topological invariant. The nonchiral QSH edge states are also different 
from the chiral edge states of the Quantum Hall insulators or Chern insulators (de- 
scribed in the previous chapter), thereby providing a new class of one-dimensional 
(ID) conductors. Interestingly the direction of the spin of the ID charge carriers is 
tied to their direction of motion. Such conductors are protected from single-particle 
backscattering (and Anderson localization) by time-reversal symmetry T. 

Unfortunately the QSH state is extremely difficult to observe in graphene due to 
the actual weakness of the spin-orbit interaction |23tl24j. In 2006, Bernevig, Hughes 
and Zhang (BHZ) predicted that CdTe/HgTe/CdTe quantum wells should host such 
a QSH state in their inverted regime [6j . Their prediction was soon followed by the 
experimental observation of conducting edge states by the group led by Laurens 
Molenkamp [71 |8]. This experimental confirmation has triggered a great deal of 
excitation in the condensed matter community [281 [29l [30l [31] . 

This chapter is organized as follows. We start by introducing the Kane- Mele 
model for graphene and the BHZ model for HgTe/CdTe heterostructures at the level 



of idealistic spin-conserving Hamiltonians (sections 3.1 and 3.2): then the QSH state 
consists in two time-reversed Chern insulators and therefore exhibits spin filtered 
counterpropagating edge modes. At this level of approximation, backscattering is 
forbidden owing to spin-conservation. Finally and most importantly, we emphasize 
that the QSH edge states are actually robust even in presence of spin-mixing terms 



owing to the protection by a bulk Z2 topological invariant (section 3.3). This in 
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variant originates from the Kramers degeneracy property applied to time-reversal 
invariant band structures of fermions. 



3.1 Kane-Mele model of graphene 

The Haldane model for spinless fermions on the honeycomb lattice (chapter [T]) and 
other Chern insulator models (chapter [2]) break time-reversal symmetry and have 
interesting topological properties akin to the Integer Quantum Hall Effect [12j . Nev- 
ertheless such Chern insulators, remain difficult to realize experimentally because 
some highly nontrivial internal magnetic background is required. In 2005, C.L. Kane 
and E.G. Mele proposed a generalization of the Haldane model that respects time- 
reversal invariance and includes the spin via the spin-orbit interaction. Their idea 
launched the field of time-reversal invariant topological insulators which has known 
a rapid expansion since then \28\ [29l [30l [3T] . 



3.1.1 Intrinsic spin-orbit coupling 

We first discuss the idealized situation of a spin-orbit coupling that still conserves 
one component of the electronic spin. In their seminal paper, C.L. Kane and E.G. 
Mele introduced the following lattice model for spinfull electrons on the honeycomb 
lattice g|: 

^ = *Y1 ^l^Cja + it2 ^ij4aiSz)af}Cj^, (3.1) 

where the Pauli matrix Sz refers to the physical spin of electrons, and the summa- 
tion over repeated spin index (a, (3) is implied. The next-nearest neighbor (NNN) 
hopping term it2VijSz describes a spin-orbit coupling between the spin direction 
= ±1 (units of h/2) and the chirality Vij of the circulating electrons. This can 
be seen as a L.S coupling where the "orbital momentum" L would be associated 
with the chirality. One can develop a very simple local picture for this spin-orbit 
coupling. 



Two copies of the Haldane insulator. Since [-ff, s^] = 0, the model Eq.(3.1) can 
be decoupled into two subsystems for spin up {sz = 1) and spin-down {sz = — 1) 
respectively. The Hamiltonian for spin- up (resp. spin-down) electrons is the Haldane 



Hamiltonian Eq.(1.25) with ip = 7r/2 (resp. = — 7r/2). Hence many properties can 
be deduced from our knowledge of the Haldane model for spinless fermions (chapter 

Insulator in the hulk. Firstly, the system is gapped in the bulk. Indeed the low- 



energy theory of the lattice Hamiltonian Eq.(3.1) is directly derived from Eq.(1.30): 



^so = AsoO-^r^s^, (3.2) 

where Ago = — 3-v/3i2- This perturbation anticommutes with the kinetic Hamilto- 
nian and therefore opens a gap at the Dirac points. For each spin specy, there is 
a mass inversion between the two valleys = — ?7i-Ka as in the Haldane model. 
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Figure 3.1: Low energy dispersion for the spinfuU Kane-Mele insulator (consisting 
in two time-reversed copies of a Haldane insulator). 



Nevertheless the global electronic system is a time-reversal invariant insulator in the 



bulk because rriKa = n^-K-a (3.1.1) 



More formally one can consider the action of the time-reversal operator T on the 



perturbation Eq.(3.2 ). For spinless electrons on the honeycomb lattice, this operator 
is T = T^K where switches the valleys and K is the complex conjugation (see 
chapter [T]). For spinfull electrons it is 

T = T^iSyK, (3.3) 

where the additional factor isy produces the reversal of the electronic spin. It is 
clear that the mass term AgoC^ r^s^ is now even under T because both and Sz 
change signs under time-reversal {txTzTx = — and SySzSy = —Sz)- 

Finally we have seen that a perturbation that respects the fundamental symme- 
tries of graphene (T and V) can open a gap owing to the presence of spin. Generally 
speaking including additional internal degrees of freedom leads to an increase of the 
number of various possible mass terms. For spinless fermions on the honeycomb 
lattice, there are only 4 possible mass terms (Semenov, Haldane and two Kekule 
distortions) and they all break some symmetry. When the spin is included, there 
are 16 different masses, some of them breaks some symmetries while others, like the 
Kane-Mele mass, respect all the symmetries [22J. 

Edge states. Each subsystem develops a spin-polarized chiral edge state whose 
circulating direction is tied to the sign of (p (see chapter [2]). Hence the global 
Kane-Mele system has a nonchiral edge state consisting in two spin filtered counter 
propagating gapless edge modes. At the level of the spin-conserving model, those 
edge states inherit the topological character of the Haldane model edge states. Of 
course, details like their dispersion relation depend on detail of the surface cut. For 
instance edge edge states along the armchair cut cross at fc = 0, whereas edge states 
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along the zigzag edge cross at k = ibvr/a (projection of the bulk Dirac points). It 
is a natural question to ask whether this topological character will survive when 
additional spin-orbit terms mixing the two copies are included. We will see below 
that the counter propagating edge states are in fact robust as long as the bulk is 
gapped and time-reversal is obeyed. 



From the toy model of the previous chapter (see section 2.3.1 ), one can deduced 
that each half-part of the QSH state will develop an edge state that circulates 
clockwise for spin up and anti clockwise for spin down. 

Spin Hall conductivity. Here we can reproduce the Laughlin argument for the 



QSH state described by Eq.(3.1). When a quantum of flux (po is added through 
the cylinder, one electron is transferred from the bottom to the top of the cylinder 
for the spin- up subsystem, leading to the e^/h charge Hall conductivity for this 
subsystem. For the other subsystem (down-spin), the one electron transfer is realized 
from the bottom to the top surface of the cylinder. As a result, the Hall charge 
conductivities cancel out (as required by T invariance), but a net spin is transferred 
from the bottom to the top of the cylinder, namely from one edge to the other. The 
corresponding spin conductance being: 

^- = 2^-U + xJ = 2^' ^^-^^ 

which is quantized. 

Hence, the intrinsic spin orbit coupling Eq.(3.2) leads to the realization of a new 
state of matter characterized by an insulating gap and a metallic edge. Nevertheless 
it is quite unlikely that the spin-orbit coupling manifests itself only through the spin- 
conserving term above. If this term is present, other terms should be present and 
mixes the spin component. Since all the previous analysis relies on the description 
of the system as two decoupled Haldane insulators, the natural question is whether 
the bulk gap, the helical edge states and the spin Hall effect will survive to the 
presence of Rashba coupling for instance. 



3.1.2 Rashba coupling 

Such a Rashba term will be present if the mirror symmetry with respect to the 
graphene layer is broken (for instance by the presence of a substrate or by a perpen- 
dicular electric field). This coupling can be described by the following low energy 
Hamiltonian [3]: 

"Hi? = \r {cTxTzSy - UySx) , (3.5) 

which mixes the two spin directions, and spoils the conservation of Sz because 
[Hr, o"z] 7^ 0. Typically spin-orbit terms are always described by some coupling 
of the momentum with the spin Pauli matrices. Here the momentum (for both in- 
trinsic and Rashba spin-orbit couplings) is implicit in the valley degree of freedom 
Tz which encodes the momentum information ibK. Note that this term is the lowest 
order (zero order) possible in the electron momentum q measured from the Dirac 
points. 
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The corresponding spectrum is [H HT] 



EaM) = a\/q2 + (Aso + I3Xr)^ + PXr, (3.6) 

where a = ±1 and /3 = ±1. As we have seen already the spectrum is gapped in 
presence of intrinsic spin-orbit couphng only. When the Rashba coupling is increased 
the gap decreases and the system become gapless when the Rashba coupling exceeds 
the intrinsic coupling. At A/j = Ago the spectrum consists in a Dirac cone and two 
gapped parabolic bands. 

The only limitation is quantitative and is related to the very weak spin-orbit 
coupling in graphene. The intrinsic SO coupling is extremely small at best 1 mK 
while Rashba can be around IK. It has been shown that SO is weak not only because 
carbon is a light element but also due to the particular arrangement of the pz orbitals 
\23\ I24|. As a consequence, in order to enhance this intrinsic spin-orbit coupling one 
should either bend the graphene layer or coat it by heavier atoms like In or . 

3.2 HgTe/CdTe heterostructures 

The Kane-Mele model is a paradigm for 2D topological insulators although the 
smallness of the intrinsic spin-orbit coupling in graphene hinders the experimen- 
tal verification. Fortunately we will see that the QSH insulator can be realized in 
materials with heavier elements and therefore larger spin-orbit coupling. Indeed 
Bernevig, Hughes and Zhang (BHZ) realized that quantum wells comprising one 
HgTe layer confined between two barriers of CdTe are good candidates for the re- 
alization of the QSH state. Those authors identified the band inversion between 
HgTe and CdTe as the crucial ingredient to realize the QSH state. Furthermore a 
minimal 4 band Dirac model (BHZ model hereafter) captures the essential physics 
of this band inversion [6J . This effective model can be derived from the 6 (or 8) band 
k.p model of 3D HgTe/CdTe complemented by an envelope method to describe the 
heterostructure. In the absence of spin mixing terms, the QSH insulator realized 
in this model also consists in two time-reversed copies of a Chern insulator. Since 
each copy has both a definite spin projection and a definite chirality, the QSH edge 
states are two spin-filtered counter-propagating modes. 

3.2.1 Band structures of HgTe and CdTe 

Three dimensional HgTe and CdTe. HgTe and CdTe are three-dimensional 
semimetal and semiconductor respectively. Both HgTe and CdTe crystallize in 
the zinc-blende band structure which is similar to diamond but with different 
atoms occupying the two inequivalent sites. HgTe is a semimetal and CdTe a 
semiconductor. For both HgTe and CdTe, the important bands are close to the 
center F of the first Brillouin zone (FBZ). Due to the strong-spin orbit coupling, 
the total angular momentum J is a good quantum number at least near k = 0. 
There is a so-called Fg band made of s-type orbital with spins up or down. There 
are p-type orbitals which are splitted into J = 1/2 (split-off F7 band) and J = 3/2 
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Figure 3.2: Edge state for a single Kramers crossing in the bulk. 



bands (rg). The Tj band can be ignored because it is far in energy and it does not 
participate in the band inversion between the band structures of HgTe and CdTe. 
At this stage, we have a 6-band model for 3D bulk semiconductors. 

Quantum wells. We now consider a thin HgTe quantum well realized between 
two identical CdTe barriers. Using the envelope function method, subband structure 
of such quantum wells have been derived from the 6-band Kane model of HgTe and 
CdTe [H]. The corresponding 6 states are denoted: \E1±) ,\H1±) and \L1±). It 
turns out that the degenerate levels electron-like \Elziz) have a band crossing with 
the degenerate heavy- hole like \Hlziz). In contrast |il±) does not participate to the 
crossing, stays far in energy, and therefore can be ignored. 



3.2.2 Effective model BHZ: massive Dirac fermion 

Continuum version. The effective model describes the low energy physics near the 
band-crossing of the 4 states \Elziz) and \Hlziz). In the absence of bulk-inversion 
asymmetry (assuming the atoms "Hg, Te and Cd " were the same atom" thereby 
transforming the zinc-blende structure in the diamond- like structure), the blocks 
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with opposite spin projections are decoupled. Therefore it is sufficient to study a 
2 block, the other block being deduced by TR symmetry. It turns out that the 
low energy dynamics of this four-band model is captured by the massive Dirac 
Hamiltonian [6] 

given in the basis order , \H1+) , \E1 — ) , \H1—)). The spin up block: 

7^(k) =eo(k)l2x2 + d(k).cr (3.8) 

is expressed in terms of the standard Pauli matrices cr acting here in the {\E1) , \H1)) 
space. The Hamiltonian ?^(k) is a Taylor expansion with respect to the in-plane 
wavevector ]i.= {kx,ky) whose coefficients are constrained by parity and time- 
reversal symmetries. 

Prom the symmetries it is possible to infer that the diagonal coupling in the 
block ^(k) are even functions of k, and the a diagonal 2x2 matrix can always 
be decomposed as el + dz(Tz- In contrast the off-diagonal terms connect opposite 
parities orbitals and should therefore be odd functions of k. Hence one has 

d{\C) = {Ak.,,Aky,M -B{kl + kl)), (3.9) 

and eo(k) = C-Dk2. 

Microscopic theory further yields he parameters A, B, C, D and M as functions 
of the quantum well geometry [B]. The parameters A ~ — 3.8eV.A and B ~ —60 
eV.A^ have a definite sign over the relevant range of thicknesses, whereas M changes 
sign at roughly d = 6.4 nm. 

The inversion between \E1) and \H1) is controlled by the sign of the mass term 
M, the QSH state being realized in the inverted regime (M < 0). In fact the 
trivial/non trivial topology of this two-band model (spin up block) is determined 
by the relative signs of B and M. Qualitatively, this can be inferred from the 
component (i^(k) = M — Sk^ which indicates that MB is positive the vector d 
points to opposite directions (north and south poles) at the F point (k=0) and 
far from it (large |k|). This configuration which looks like a meron in k— space is 
nontrivial. In contrast if MB < 0, then the texture d(k) has no winding. Since 
material parameters lead to 5 < 0, the nontrivial regime corresponds to the inverted 
regime M < 0. 

The chemical potential C determines the electronic filling of the bands which 
can be electrostatically tuned by the action of a distant metallic gate. 

Regularized lattice version. The above model is already simplified and has the 
advantage of being quantitative (values of parameters derived from a more micro- 
scopic theory). Nevertheless being ak.p expansion it is only valid locally near k = 0. 
Of course one cannot solve the full band structure in the whole FEZ, but one can 
imagine a simple idealized model which is defined over the whole FEZ and has the 
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same low k expansion as effective model Eq. (3.9), namely: 



d(k) = {Asmkx,Asmky,M - 2B{2 - cosk^ - cosky)). (3.10) 

This can be represented as a square lattice with four states at each site: ■01,3 = s is 
s-type orbital with up/down spin, while V2,4 = Px i ipy-spin orbit coupled orbital 
with up/down spin. 



3.2.3 Edge states of the BHZ model: interface with the vacuum 

We now investigate the edge states of HgTe/CdTe quantum wells using the BHZ 
model which neglect the spin-mixing BIA and SIA terms. Then it is sufficient to 



consider the spin-up block 'H(k) in Eq.(3.7). For simplification, we neglect the eold 



term which is not relevant for the topological properties. 

We consider an interface along y— axis between the HgTe/CdTe well (x > 0) and 
vacuum (x < 0). Owing to translational invariance along y— axis, the momentum 
ky is a good quantum number whereas kx should be replaced by the real-space 
derivative dx- In the half-plane x > 0, the wave function then obeys a sec ond-order 
differential equation: 

{-iAaxdx + Akyay + {M + Bdl - Bk'^)a,) ^ = E^, (3.11) 

with the strict boundary condition ^{x = 0) = 0. We search for a zero energy mode 



at ky = 0. Hence we start by solving Eq.(3.11) with ky = and E = 0: 

{-iAaxdx + {M + Bdl)a,)^ = 0, (3.12) 
which can also be written as: 

- Bdl"^ - Auydx"^ = M"^. (3.13) 
If we inject an exponential solution ^' = ipe^^, one finds 

(M + 5A2)$^ = -AayX<^^, (3.14) 
and therefore the solution are eigenmodes of ay with decay lengths given by: 

X:^ = ^ , (3.15 

^ 2B ' ^ ^ 

where in this equation ay has to be understood as an index cjy = it referring to the 
eigenvalue of the matrix ay. Then if we solve the problem in the half-plane x > 0, 
one should pick up the solutions with ReA < which have physical asymptotic 
behavior at x — )• oo. Then two scenarios are possible. If two solutions have the 
same Uy-polarization, then it is possible to satisfy the condition ^(x = 0) =0 and 
the solution reads: 

^{x) = ^^(e^^--e^n, (3.16) 
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where ay is the chosen polarization. Since A^^.A*^" = M/B, this situation is reahzed 
if M/B > 0. In contrast, when M/B < 0, the two sohitions having the right 
asymptotic behavior belongs to orthogonal polarizations, and it is impossible to 
cancel the wave function at a; = by the superposition of two orthogonal spinors. 

In terms of the parameters, the condition for having an edge state turns out to 
be > 4MB > 0. The signs ^4 > and B < are fixed over a broad range of 
HgTe width d, whereas the "mass" M changes it sign as a function at dc = 64 nm. 
According to the BHZ model, there is a zero mode E = at ky = when M < 0. 

Dispersion and Hamiltonian of the edge state. For M < 0, we can restore a 
finite ky and notice that the zero energy solution is still an eigensolution at energy 
E = ayAky if one neglects the terms Bk^Uz (which is always possible at very low 
ky). The (7j^ = ± indicates the polarization of such an edge mode. Note that such 
an edge mode is chiral which was expected because it is the edge mode of the Chern 
insulator described by the spin-up block of the BHZ model. The full model has an 
edge mode carrying opposite spin, opposite Gy polarization and circulating in the 
opposite direction. 

3.2.4 Interface between two BHZ insulators with different masses 

Another model we may want to solve is the interface between two BHZ insulators 
characterized by distinct values of the mass. The corresponding equation being: 

{-iAaA + AkyGy + (M(x) + Bdl - Bkl)(Jz) * = E^, (3.17) 

with M{x) = MiQ{-x) + M2Q{x). 

Now we can try to construct the edge states solutions for a semi-infinite HgTe 
well extending in the half-plane ?/ > 0. Boundness requires to choose the two 
exponential with negative A while the boundary condition is ^{x,y = 0) = 0. Note 
that in the particular model of HgTe wells near the transition point, parameters A 
and B are both negative. Thus Ai is always negative whereas the sign of A2 is the 
sign of —MB. We now discuss the cases MB > and MB < 0. 

Conclusion: We find no edge state betw^een two nontrivial (or two 
trivial) insulators. We find a gapless edge state between a trivial and a 
nontrivial insulators. Note that the edge state consists in a single expo- 
nential on the trivial side while it is characterized by two exponentials 
on the nontrivial side. 

3.3 Helical edge states 

In the previous examples of spin conserving models, we have seen that QSH insu- 
lators consist in two copies of a Chern insulator (or QAH state), each copy being 
associated with a spin orientation. Since the QAH has (minimaly) a single chiral 
edge state, we can deduce that the QSH state will have two spin-filtered counter 
propagating states. Due to spin-conservation, there is no backscattering coupling 
those states. Nevertheless one may question if the edge states pertain for more real- 
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istic models that include the spin-mixing unavoidably present in any material with 
stron spin-orbit coupling. 

The crucial question is the robustness of this helical state in presence of spin mix- 
ing. Numerical indications that the edge states are preserved as long as the bulk is 
gapped. This paragraph tries to explain the general reason for this protection which 
is Kramers degeneracy. We aim at describing the implications of Kramers degener- 
acy for Bloch electronic band structures, both at the level of eventual edge states 
and at the bulk level. Finally the name quantum spin insulator can be mislead- 
ing. This is a new insulating state that is different from the doped semiconductors 
exhibiting QSH effect (metals versus insulators). Nevertheless if the Fermi level is 
raised into the 2D conduction or valence bands, we obtain such states. 



3.3.1 Kramers degeneracy 

The time-reversal operation T is anti-unitary. Then = ±1 depending on the spin 
of the system. We are interested in fermionic systems, and then = — 1. Then 
one can prove that for a time-reversal invariant system, characterized by [if, T] = 0, 
any state is always degenerated. Indeed if 99 is a stationary state with energy E, 
namely Hip = Eip, then T^p is also an eigenstate of H with the same eigenvalue. 
Moreover Tip is a different quantum state than ip precisely because of = — 1. 
Indeed if T(p and ip were the same state, then one could find a complex number c 
realizing Tip = ap and therefore T^ip = c*Tip = \c\'^<P- The last equality leads to 
the contradiction IcP = —1 because = —1. 



3.3.2 Configuration of edge states 

The spirit of this section is the following. Let us assume the existence of some ID 
edge states within the bulk gap. Then we shall examine the implications of time- 
reversal symmetry on the dispersion of such edge states The dispersion is 
defined with respect to a momentum parallel to the border of the insulator, denoted 
/c, which lives in the segment [— vr/a, vr/a]. The time-reversal invariant points are 
k = Q and k = Ti/a (which is the same as /c = —vr/a). 

A single Kramers doublet. The minimal situation corresponds to 1 single 
Kramers doublet forming at the TRI points that we can take at = for simplicity 



as in Fig. 3.3 When k is varied from to vr/a say, the degeneracy is lifted and 
the two levels can move inside the bulk gap. Assuming a single Kramers doublet 
for this surface state, the splitted Kramers partners have to merge into the bulk 
bands either the conduction or valence bands. Then we have two possible scenarios 
depending whether the Kramers partners merge two the same bulk band or to 
distinct bands. 

Two Kramers doublets. We now imagine a slightly more complex situation with 
Kramers doublets at fc = and k = ir/a within the bulk gap. We also assume 
that the states belonging to those two doublets will become intricate in some way 
(otherwise we are left with the previous situation of two independent decoupled 



Kramers pairs). Then there are two possible scenarios illustrated in Fig. 3.3.2 
First, the two partners can stay in the gap region and merge again together at 
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Figure 3.3: Edge state for a single Kramers crossing in the bulk. 



k = 'K/a (Fig. 3. 3. 2, a). Second, one of the partner can merge into the bulk valence 



band while the other will stay away from the band edge and will be joined by a state 
arising from the conduction band (3.3.2 b). In the first case, the Fermi level in the 
gap will cross either two or zero pair(s) of edge states. This is similar to a standard 
ID conductor. Moreover if the Kramers degeneracy at TRI points are moved both 
up to conduction band (or both down to valence band) the edge states disappear. 

In the second case, the FL intersects only a single pair. The two counter propa- 
gating modes at the Fermi level are each non degenerate and carry a given Kramers 
parity. This is half of a usual ID system. Moreover if we move the two Kramers 
doublet to the bulk bands (in any way) , there will always remain a pair of Kramers 
partners crossing the bulk gap. 



3.3.3 Absence of single-particle backscattering 

We have seen that the crossing of the time-reversed states is protected by T at the 
time-reversal invariant points of momentum space. Now we investigate the relation 
between time-reversed states, typically labelled by k and — k, at any momentum k. 
It can be proved that single-particle processes between such states are forbidden. 
This extends dramatically the absence of backscattering discussed in the previous 
sections on spin-conserving models. 

This can be also proved in a very general way general way. Let us consider two 
time-reversed states (p and ijj = Tip. Then we can show that the matrix element 
{tplHip) is always zero provided the system is time-reversal invariant [H, T] = and 
fermionic = — 1: 
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a a a a 

Figure 3.4: Edge state for two Kramers crossing in the bulk. 



{iPlHip) = {TiplHip) = {TH^\T\) (3.18) 

= -{THip\ip) using = -1 (3.19) 

= -{HT(p\<f) using from[r, H] = (3.20) 

= -{T^\H\^) (3.21) 

There is another proof using second-quantization and based on the transformation 
rule of creation/annihilation operators for one- half particles (29j : 

TcktT-^ = (3.22) 
Tc-hT'^ = -Ckt (3.23) 

Due to those properties the Hamiltonian describing single-particle backscattering 
for a single pair of Kramers partners [29] : 

^back = 4.c-k-s + H.c, (3.24) 

k,s 

is odd under T: 

T-ffbackT"^ = — -ffback- (3.25) 

Finally, we note that two-particle processes are allowed. Recently it was dis- 
cussed that they can be generated microscopically by the interplay of Rashba cou- 
pling and electron-electron interactions. 



36 



3.3.4 Z2 topological invariant 

The central idea is to implement the condition of time-reversal invariance in the con- 
text of systems where spin-orbit induces a band inversion and a related topological 
state. If the edge states are present, then TRI enforces that their dispersion curves 
should meet at the TRI points of the projected BZ. The bulk states should also be 
constrained by time-reversal symmetry. When TR is obeyed the TKKN invariant is 
zero. Kane and Mele shown that there is an other invariant which can be = 0, 1 
only. There is a formula for this Z2 topological invariant which can be expressed as 
an integral over the whole bulk BZ. For systems having a center of inversion, the 
additional V symmetry allows to simplify this formula and to calculate the topolog- 
ical invariant z/ = has a simple product of the parities of occupied states at all the 
inequivalent TRI points of the FBZ. 

In the QH effect (and in Chern insulators) there is a correspondence between 
the value of the bulk topological invariant of the occupied bulk bands (the TKKN 
or Chern number) and the number of chiral edge states. By qualitative consider- 
ations we have found that there should be a Z2 classification of the edge states of 
time-reversal invariant insulators related to the even/odd effect in the number of 
edge modes crossing the Fermi level. In presence of inversion symmetry, there is 
a simplified formula to evaluate the Z2 invariant v and determine whether a given 
band structure is a nontrivial topological insulator or a trivial band insulator. The 
formula is simply [49j: 

(-i)'^=n n ^2™, (3.26) 

kinv m=\,..,N 

where (^2ra = C2m+i is the parity of the Bloch wave function at k = kinv in the 
occupied band labelled by 2m. 

Conclusion 

The Quantum Spin Hall (QSH) effect is a property of certain two-dimensional elec- 
tron systems with spin-orbit coupling: the bulk of the system is electrically insulat- 
ing, while a conducting "helical metal" exists at the boundary in which electrons of 
opposite spins move in opposite directions [HIS]. The prediction [6j and the observa- 
tion (TIE] of the Quantum Spin Hall (QSH) state in mercury telluride (HgTe/CdTe) 
heterostructures have triggered a great deal of excitation in the condensed mat- 
ter community j 28| [29l |3T] since the QSH state realizes a two dimensional (2D) 
topologically ordered phase in the absence of magnetic field. 

Finally we mention that the QAH and QSH phases might be generated sponta- 
neously by some types of local Hubbard correlations [50^ l5T] even in material with 
light elements with weak spin-orbit couplings in the absence of interactions. 
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